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The  objective  of  this  study  is  to  determine  the  effect  of  clouds  on  clear  Landsat  pixels  that 
are  adjacent  to  cloudy  pixels.  A  partly  cloudy  scene  over  a  well-defined  ocean  background 
was  analyzed  and  the  effect  of  clouds  on  the  clear  pixels  as  a  function  of  cloud  amount  was 
determined.  The  amount  of  scattering  that  occurs  is  proportional  to  the  amount  of  cloud 
'  edges  in  the  vicinity,  which  is  proportional  to  cloud  amount.  For  each  non-cloudy  pixel,  the 
'  cloud  fraction  was  computed  in  a  7x7  box  using  a  simple  threshold  technique,  and  the 
:  average  brightness  in  that  box  was  compared  with  the  brightness  you  would  expect  to 
measure  if  there  was  no  scattering  taking  place: 

Obs  =  (l-n)Ocean  +  (n)Cloud 

where  Obs  =  observed  Landsat  grayshade,  Ocean  =  dear-column  ocean  brightness 
value.  Cloud  =  grayshade  value  of  a  Landsat  pixel  that  Is  completely  filled  with  a  cumulus 
cloud,  and  n  Is  the  cloud  amount  in  the  7x7  box.  Any  difference  from  this  behavior  is  due 


to  scattering. 

Adjacency  effect,  Cloud  fraction.  Cloud  scattering 
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1.  INTRODUCTION 

The  objective  of  this  study  is  to  determine  the  effect  of  clouds  on  clear 
LANDSAT  MSS  pixels  that  are  adjacent  to  cloudy  pixels.  Water  droplet 
clouds  are  efficient  scotterers  of  visible  and  near-infrared  solar  radiation. 
Scattering  of  radiation  by  clouds  -within  a  particular  field  of  -view  can 
affect  nearby  pixels  that  do  net  contain  clouds.  A  partly  cloudy  scene 
over  a  well-defined,  uniformly  dork  ocean  background  -was  anedyzed 
and  the  effect  of  douds  on  the  adjacent  clear  pixels  was  determined  os 
a  function  of  cloud  amount. 

2.  ANALYSIS  METHOD 

Londsat  MSS  data  from  a  scene  in  the  Bering  Sea  on  26  August  1980 
•was  chosen  for  this  analysis.  The  scene  contains  a  group  of  small-scale 
cumulus  clouds.  Solar  elevation  is  sufficiently  high  for  good  scene 
illumination,  and  the  uniformly  dark  Bering  Sea  serves  as  a  useful, 
•wen-behaved  background  o-ver  which  to  investigerte  the  attributes  of 
scattering  by  cloud  edges  into  adjacent  fields  of  -view.  The  over-'water 
background  also  simplifies  the  complex  effects  of  doud  shadows.  Qoud 
imagery  for  the  selected  scene  is  shown  in  Figure  1. 

The  magnitude  of  scattering  by  douds,  called  file  "doud  adjacency 
effect,"  is  proportional  to  the  amount  of  doud  surface  area  exposed  to 
incident  solar  radiation  in  the  -vicinity  of  a  field  of  -view.  Although  not 
precise,  doud  fraction  is  used  in  this  study  to  specify  the  amount  of  doud 
surface  area  near  a  given  pixel.  Cloud  fraction  is  determined 
simultaneously  for  each  MSS  image  band  4-7  by  applying  a  threshold 
brightness  below  which  no  douds  appear  to  be  present  in  the  image. 
MSS  Bond  7  derta  are  useful  for  determining  this  threshold  since 
scerttering  by  the  atmosphere  and  by  douds  is  minimal  at  the  longer 
near-infrared  wavelengths.  The  clo-ud  mask  was  generated  for  each 
pixel  and  for  each  MSS  band  in  the  scene  by  flagging  all  pixels  -with  a 
bond  7  grayshade  greater  than  15.  Figure  2  contains  this  doud  Tnn<;k: 
for  the  images  in  Figure  1. 
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Cloud  fraction  “in  the  vicinity"  at  each  pixel  is  then  computed  by 
counting  the  number  of  cloudy  pixels  in  the  sunounding  7x7  box,  and 
dividing  through  by  the  total  number  ot  pixels.  Qoud  fraction  is  used  os 
a  measure  of  the  amount  of  cloud  surface  area  that  is  exposed  to 
incident  solar  radiance.  Hence,  fraction  is  a  measure  of  the  degree  to 
which  cloud  scattering  takes  place. 

If  the  cumuliis  clouds  from  one  pixel  have  no  effect  on  the  upwelling 
radiance  measured  in  an  adjacent  pixel,  then  the  dependence  of 
radiance  on  doud  fraction  (cloud  surface  area)  in  a  7x7  box  is  linear, 
and  would  be  described  as 

^OBS  "  ^^CLD  ^CLR  “  ^^CID~^CLR^^  ^  ^CLE  ' 

where  Iq^j,  p  ore  the  observed  Landsat  MSS  radiances  for  a 

jxirtly  doudy,  completely  cloud-filled,  and  doud-free  (water)  pixel, 
respectively,  and  where  n  is  doud  fraction  (0<n<l).  Any  deviation  from 
this  depjendence  represents  the  effed  of  scattering  into  the  field  of  view 
by  doud  edges  in  adjacent  pixels. 

A  plot  of  radiance  (which  for  Landsat  MSS  bonds  4-7  are  linearly 
proportional  to  grayshade)  as  a  function  of  doud  amount  car  be 
produced  from  the  image  data  for  each  MSS  band.  For  each  pixel  in  the 
scene  the  local  doud  fraction  n  is  computed  for  the  surrounding  7x7  box. 
Then  the  grayshade  for  that  pixel  is  tcbuloted  as  a  function  of  doud 
fraction.  Once  completed  for  the  entire  image,  the  overage  grayshade 
for  each  doud  fraction  is  then  computed  and  compared  to  the  linear 
dicqae  that  is  assumed  by  Eq.  (1).  This  plot  is  shown  in  Figure  3. 

3.  RESULTS  AND  SUMMARY  REMARKS 

fri  Figure  3  it  is  seen  that  the  doud  adjacency  effect  is  both  nonlinear 
and  wavelength-dependent  for  MSS  bands  4  and  5.  For  bands  6  and  7 
the  effect  is  essenticdly  non-existent;  that  is,  radiance  is  linear  with  doud 
fraction.  MSS  bands  at  the  lower  wavelengths  are  most  strongly 
affected  by  scattering.  Also,  the  adjacency  effed  is  strongest  for  n  =  0.7. 
The  vertiodl  difference  between  the  dotted  straight  line  in  Figure  3  [the 
linear  assumption  as  defined  by  Eq.  (1)]  and  the  curves  for  bonds  4  and 
5  quantify  the  effed  of  doud  scattering  into  non-cloudy,  adjacent  fidds 
of  view. 

Curves  such  os  Figure  3  may  prove  useful  in  removing  the  scattering 
effects  of  clouds  from  Landsat  data  in  order  to  more  accurately  retrieve 
the  refledance  attributes,  such  as  vegetation  indices,  of  doud-free 
surfaces.  Also,  if  doud  fraction  is  the  retrieved  parameter  of  interest, 
adjocerrey  effects  when  not  propjerly  accounted  for  generally  lead  to  an 
overestimation  of  doud  fraction. 
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Future  work  wifl  include  reconstnrcting  the  Bering  Sea  scene  using 
the  results  plotted  in  Figure  3  to  determine  whether  the  corrected 
doud-free  pixels  have  a  unitoim  dork  appearance.  More  challenging 
studies  include  analysis  of  scenes  with  non-water  backgrounds  where 
shadowing  effects  significantly  complicate  the  satellite-measured 
reflectances. 
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Figure  1.  Londsat  MSS  Bond  4  (0.5-0.6  |im.  Upper  Left),  Bond  5 
(0.6-0.7  pm.  Upper  Right),  Bond  6  (0.7-0.8  nm.  Lower  Left),  and  Bond  7 
(0.8-1. 1  pm.  Lower  Right)  Images  for  the  Bering  Sea,  26  August  1980. 
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Figure  2.  Landscrt  MSS  Images  of  Figure  1,  •with  the  Qoud  Mask 
Applied.  Bright  White  Pixels  ore  Qoud-Filled;  cdl  Others  ore  Qoud-Free. 


Cloud  Fraction  (Z) 

Figure  3.  Londsat  MSS  Bonds  4  and  5  Normalized  Radiances  as  a 
Function  of  Cloud  Fraction.  The  Dotted  Straight  Line  Denotes  the 
Non-Scattering  Condition  Specified  by  Eq.  (1). 


